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Abstract 



For many years, ^He typified Bose-Einstein superfiuids, but recent ad- 
vances in dilute ultra-cold alkali-metal gases have provided new neutral su- 
^ \ perfiuids that are particularly tractable because the system is dilute. This 

O ' chapter starts with a brief review of the physics of superfiuid ^He, followed by 

the basic ideas of Bose-Einstein condensation (BEC), first for an ideal Bose 
gas and then considering the effect of interparticle interactions, including 
^ \ time-dependent phenomena. Extensions to more exotic condensates include 

^^ \ magnetic dipolar gases, mixtures of two components, and spinor condensates 

that require a focused infrared laser for trapping of all the various hyperfine 
cn ■ magnetic states in a particular hyperfine F manifold of mp states. With 

cn . an apphed rotation, the trapped BECs nucleate quantized vortices. Recent 

i:^ \ theory and experiment have shown that laser coupling fields can mimic the 

effect of rotation. The resulting synthetic gauge fields have produced vortices 
in a nonrotating condensate. 
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1. Introduction to Bose superfiuids: brief review of superfiuid '*He 

Superfiuid ^He has played a central role in understanding the physics of 
Bose-Einstein condensation. It has the unusual property of remaining liquid 
to T = K under atmospheric pressure. This behavior reflects its weak 
interatomic potential and the relatively large zero-point energy (because of 
the small mass). The common helium isotope '^He has two electrons, two 
protons and two neutrons and thus acts like a boson under exchange of two 
such atoms. Below a characteristic temperature Tx ~ 2.17 K, it undergoes 
a phase transition from a normal weakly viscous fluid to a superfiuid with 
remarkable properties (for a review of this fascinating system, see |l|). 

1.1. two- fluid model 

Superfiuid ^He has remarkable hydrodynamic properties. It can flow 
through fine channels with no pressure drop, which suggests that it has zero 
viscosity, yet a direct measurement of the viscosity through torque on a ro- 
tating cylinder yields a value comparable with that for the normal phase 
above Tx. This seemingly contradictory behavior was explained phenomeno- 
logically by a two-fluid model, in which the low-temperature phase has two 
interpenetrating components: (1) a superfiuid with zero viscosity and irrota- 
tional flow velocity Vg with V X v^ = and (2) a normal (viscous) fluid with 
general velocity t;„. Each component has its own temperature-dependent 
mass density ps and p„, with ps + pn = P (the total mass density). The 



superfluid component has p^ = at Tx and p^ = p at T = K. The normal 
density p„ has been measured by the viscous drag on a set of oscillating disks 
that respond only to the normal component; the superfluid density is then 
given by ps = p - p„. 

Superfluid ^He can exhibit persistent currents in a multiply connected 
geometry like a torus or in a porous random medium such as a compacted 
powder. If the container and fluid rotate above Tx and are then cooled 
below Tx the superfluid will continue to rotate at essentially the original 
angular velocity Q. If the container is slowly brought to rest, in many cases, 
the superfluid continues to rotate (in contrast to the behavior of a classical 
viscous fluid that would eventually come to rest). Such persistent currents 
typify both superfluid ^He and superconductors (where the charge makes 
the detection of such electrical currents relatively easy). There is a critical 
velocity Vc for superfluid flow that is set by the energy of the lowest excited 
state; Vc gives an upper limit to the magnitude of persistent currents in 
confined geometries. For a channel of lateral dimension ~ d, the typical 
observed critical velocity is of order Vc ~ h/Md, where M is the atomic 
mass. 

1.2. quantized circulation 

For low flow velocity, liquid ^He is effectively incompressible with ^ -Vg = 
0. Since Vg is also irrotational, it can be described with a velocity potential 
$ that satisfles Laplace's equation, with Vg = V$. For a quantum fluid, 
we can be more speciflc, with $ = {h/M)S, where M is the atomic mass 
and S is the phase of a macroscopic one-body wave function characterizing 
the quantum-mechanical condensed state of superfluid ^He at T = K. The 
irrotational condition on Vs apparently suggests that the superfluid would 
not rotate at T = K, but experiments showed that the rotating superfluid 
had a parabolic meniscus independent of temperature (like that of a classical 
viscous fluid). In 1949, Onsager proposed without explanation that the su- 
perfluid circulation k = j dl-Vg was quantized in units of h/M. Subsequently 
Feynman [^ suggested that rotating superfluid '^He has a uniform array of 
quantized vortex lines with singular vorticity n at the center of each lineo 
The velocity potential for a single vortex line is proportional to the azimuthal 
angle 0, which ensures that the quantum-mechanical wave function is single- 

^Donnelly |3[ gives good account of this intriguing early history in Sec. 2.3. 



valued and reproduces Onsager's quantized circulation. Feynman also chose 
the areal vortex density n^ to mimic solid-body rotation v^h = fi x r, yielding 
n^ = MQ/ttH, where Q is the external angular velocity. 

1.3. quasiparticles and the Landau critical velocity 

Landau [4, |5|] introduced a quasiparticle model to describe the normal 
component (for more accessible accounts, see [ll, l6|). Specifically the ele- 
mentary excitations act like phonons at long wavelengths with a dispersion 
relation ujk = sk, where s ~ 240 m/s is the speed of compressional sound. 
For shorter wavelengths comparable with the interparticle spacing in liquid 
^He, the dispersion relation has a local minimum associated with what are 
called rotons 

Uk ^ h n- —, (1) 

" h 2M* ' ^ ^ 

where A/ke ^ 8.7 K is the roton gap, ko ^ 1.9 x 10^° m"^ and M* ^ 0.16M 
for ^He. This model of independent phonons and rotons allows a straight- 
forward calculation of the temperature-dependent normal- fiuid density p„. 
More importantly, it also predicts a critical velocity as follows. Imagine a 
macroscopic object moving through the superfiuid at T = K. Conservation 
of energy and momentum indicates that this object cannot lose energy by 
creating a quasiparticle unless its speed exceeds a critical velocity Vc given 
by the minimum value of Uk/k: 



k 



(2) 
For superfiuid "^He, rotons determine this value as Vc ~ A/hko ~ 50 m/s. 



2. Three-dimensional ideal Bose gas 

This section provides a brief review of the ideal Bose gas, both the familiar 
uniform case and for a gas confined in a harmonic trap as in most experiments 
with ultracold atomic vapors (good general references are [7|, |8|, |9|). 

2.1. qualitative picture of Bose-Einstein condensation in a uniform system 
Consider a uniform ideal gas with number density n and volume per 
particle n~^. One important characteristic length is the interparticle spacing 
~ n~^^^. Suppose that the system is in thermal equilibrium at temperature 
T. If M is the particle mass, the mean thermal momentum is pt ~ ^/MkBT. 



The de Broglie relation then gives the thermal wavelength Ay ~ h/^MkBT, 
which provides another characteristic length for this uniform ideal gas. 

In the classical limit /i — )■ or at high temperature, the thermal wave- 
length becomes small compared to the interparticle spacing n~^/^, which 
means that the effect of quantum diffraction and interference is negligible. 
This limit is equivalent to ray optics for light when the wavelength is very 
small compared to the typical dimension. 

As T falls, however, the thermal wavelength eventually becomes compara- 
ble with n"^/^, and quantum interference between nearby particles becomes 
crucial. It is convenient to introduce the dimensionless phase-space density 
n\\, which is small in the classical limit. When nA|i is of order unity, how- 
ever, the ideal gas changes character. An ideal gas of fermions undergoes a 
crossover from the classical ideal gas to a degenerate Fermi gas at the Fermi 
temperature Tp given by the qualitative relation n\\^ ~ 1. Equivalently, 
ksTp ~ f?r?''^ jM^ is the usual expression in condensed-matter physics. For 
n ~ 10^^ m^^ and M equal to the electron mass, the Fermi temperature 
Tp for electrons in a metal is of order 10^ K, whereas liquid ^He has similar 
number density but particle mass ~ 10'* times larger, yielding Tp ~ 1 K. 

In contrast, an ideal gas of bosons in three dimensions undergoes a sharp 
phase transition at a temperature T^ given by the similar criterion: n\^ ~ 
1, or equivalently ksTc ~ W"n^l^ jM^ similar to that for a Fermi gas. For 
liquid ^He, the number density and mass are comparable to those for ^He, 
so that the transition temperature is of order 1 K. Nevertheless, the physics 
of the Bose-Einstein condensation is quite different from the onset of Fermi 
degeneracy, and we now explore the specific cases of a Rose gas in a box and 
in a general harmonic trap. 

2.2. quantitative description of ideal Bose gas in an external potential Vtr 

Consider an external trap with potential Vtr, and an associated complete 
set of quantum- mechanical single-particle states with energies ej. Assume an 
ideal Bose gas in equilibrium at temperature T and chemical potential /x in 
this external potential. The mean occupation number of the state j is 

where (3~^ = ksT, and /(e) = {exp[/3(e— /x)] — 1}~^ is the usual Bose-Einstein 
distribution function. The mean total particle number and mean total energy 



are given by 

N{T,f,) = J2fi^,), (4) 

j 

S(T,/i) = ^6,/(e,). (5) 

j 

Formally, one can invert Eq. (jlj) to express the chemical potential as /i(T, A^), 
thus giving the energy in terms of the more familiar variables T, N. 

Since we shall deal with many states, it is convenient to introduce the 
density of states g{e) = ^ • 6{e — e^) where the sum is over all distinct single- 
particle states (some of which may be degenerate). This applies for any 
particular trap potential. Equation (^ can then be rewritten as 

N{T,i^) = jdegie)fie). (6) 

In the classical limit for fixed A^, the chemical potential is large and nega- 
tive [since rij <^ 1 in Eq. ([3])], but as the temperature decreases, the chemical 
potential increases and eventually approaches the lowest single-particle en- 
ergy eo. This equality fi(Tc, N) = eo defines the critical temperature for the 
onset of Bose-Einstein condensation (BEC). For T < Tc, the chemical poten- 
tial cannot increase, because the distribution function /(eo) would become 
singular. In contrast, Eq. (E]) continues to decrease for T < Tc. The expla- 
nation is that a macroscopic number of particles starts to occupy the lowest 
single-particle state with occupation number No{T). The temperature de- 
pendence of A'"o(T) follows from conservation of particles A^ = Nq{T)+N'{T), 
where 

defines the total number of particles not in the condensate. Note that this 
integral is finite only if g{eo) vanishes; otherwise N'{T) diverges and BEC 
cannot occur. 

2.3. ideal Bose gas in three-dimensional box with periodic boundary condi- 
tions 

Imagine a cubical box of length L on a side. Plane waves have the familiar 
normalized wave function tpir) = V~^l'^ exp{ik ■ r) with V = L^ and single- 
particle energy e^ = h^k"^ /2M. The states obey periodic boundary conditions 



for the specific wave vectors k = {2ti/L) {nx,ny,nz), where Uj is any integer. 
The lowest single-particle state is uniform with k = and cq = 0. The 
corresponding density of states is readily found to be 

which vanishes at e = 0. The condition for the onset of BEC yields the 
familiar condition for the phase-space density n\\^ = C(|) ~ 2.612 at Tc. 

Below Tc, the fraction of particles in the excited states (those with finite 
momentum) is given by 



N \T, 

Conservation of total particles then yields the fraction of condensed particles 
in the zero-momentum single-particle ground state 

'm^i-dY', (10) 



N " \fj 

which increases continuously from at Tc to 1 at T = K. 

2.4- ideal Bose gas in three-dimensional harmonic trap 

Typical experiments on dilute alkali-metal gases rely on magnetic traps 
that vary quadratically with the distance from the origin (see Chap. 1). 
Hence we consider a general three-dimensional harmonic trap potential 

VUr) = \M {ujIx^ + a;y + ulz^) (11) 

with the familiar single-particle energies 



e 



f^x'^y'^z 



h {n^uj^ + riyUJy + n^Uz) + eo, (12) 



where nx,ny,nz are non-negative integers and eo = ^h{ux + Uy + Uz) is the 
zero-point energy in this harmonic trap. The ground-state wave function 
ilJo{x,y,z) is a product of three Gaussians with characteristic dimensions 
dj = ^Jh/MuJj, where j = x, y, z. If the sums in the density of states are 
approximated by integrals (which holds for e ^ eo), the resulting density of 
states becomes 



where Uq = (w^cuj^w^)^/^ is the geometric-mean frequency. 

The onset of BEC in this harmonic trap occurs for fi{Tc, N) = eo, with 
the transition temperature 

ksTc = [C{3)]-^/^huoN^^^ ^ 0.94 huj^N^'^, (14) 

where ^(3) ~ 1.212. Below Tc-, the number of particles in the excited states 
decreases like N'{T)/N = {T/Tc)^. The remaining particles occupy the 
ground state so 

'^"'^'-l-f^V. (15) 



The condensate fraction grows from at Tc to 1 at T = K as in the case of a 
uniform Bose gas, but with a temperature dependence different from that in 
Eq. ([9]). The single-particle ground-state wave function is that of a harmonic 
oscillator, ipQ{x,y,z) oc J^ •exp(— a;^/2(i^). Note that the nontrivial trapping 

potential VJ;r introduces a new characteristic oscillator length do = y^h/Mcuo 
that does not appear for the uniform ideal Bose gas. For a noninteracting 
gas (or small numbers of atoms) the zero-temperature condensate density 
has a Gaussian profile nQ{r) = N\ipo{r)\'^. For an interacting gas the size of 
the condensate is larger than the harmonic oscillator length as shown in Sec. 
14.11 but Eq. f lTSj) remains a good approximation since it is the Bose-Einstein 
statistics that determine T^, not the interactions. 

Images of the temperature-dependent condensate \1(3\ provided clear ev- 
idence for the formation of a BEC in dilute cold alkali-metal gases (see, for 
example, [7|, |8|, |9| ) . Typical atomic traps have t/o ~ a few /im, and A^ ~ 10^, 
leading the transition temperatures of order 100-1000 nK, depending on the 
atomic mass. 

2.5. ideal Bose gas in harmonic traps with reduced dimensions 

The shape of the confining potentials for ultracold atoms (created by ap- 
plied magnetic fields, laser light or other external fields) can be engineered 
to be highly asymmetric so that certain degrees of freedom are frozen out. 
Strong confinement in two directions, u^ <^ cOx — ^y, creates a long thin tube 
of atoms and such one-dimensional systems have some fascinating properties. 
For example in a Tonks-Giradeau gas (see Chap. 4) bosonic atoms act like 
fermions (in some but not all ways). This behavior arises when the confine- 
ment is so tight that atoms cannot pass each other, and this one-dimensional 
system of impenetrable bosons resembles beads on a string, or cars in a traffic 



jam — atoms that cannot be at the same position along the system look as if 
they obey the Pauli exclusion principle. 

Conversely for very tight confinement in one direction, Uz ^ cOx — ^y in 
Eq. flTTj) . the cloud of atoms has a oblate (pancake) shape; when the temper- 
ature is low enough k-gT <^ hwz the atoms move only in the radial plane giv- 
ing a two-dimensional gas, e.g. for lOzj'^'^ = 2 kHz this requires T ^ 100 nK 
which is readily achievable. Two-dimensional behavior has been investigated 
in films of superfluid helium but an important difference is that the helium 
system corresponds to a square well (with a flat potential between the walls) 
whereas in ultracold gases the atoms are usually in a harmonic potential — 
these two situations have different densities of the states [as we have seen 
for the three-dimensional cases in Eq. ([H]) and Eq. 013p ]. A two-dimensional 
gas in a uniform potential is a special case: the density of states is indepen- 
dent of energy so the integral that gives N' diverges in Eq. ([7]) — this means 
that BEC does not occur in this case. Nevertheless interesting phenomena 
occur in a two-dimensional gas in the quantum degenerate regime such as 
the Berezinskii-Kosterlitz-Thouless transition, when the density of atoms per 
unit area is sufficiently low, and the physics related to this is discussed in 
Chap. 4. Other trapping geometries have also been demonstrated, one of the 
most interesting being the creation of ring-shaped clouds where there can be 
persistent flow around the loop, as for superfluid helium in similar multiply 
connected systems. Another case is that of rapid rotation about the axis of a 
cylindrically symmetrical potential which is described described in more de- 
tail in Sec. |9l centripetal acceleration causes the gas cloud to bulge outward 
along its equator so that it can become flattened into the 2-D regime, i.e. 
effectively the tightness of trap is reduced in radial direction in the rotating 
frame, as shown in Eq. fIBU]) . 

3. Energies and length scales for interacting Bose gas 

The inclusion of a two-body interaction potential U{r — r') leads to sev- 
eral new features that play a crucial role in the description of realistic dilute 
cold Bose gases. Typically, the interactions have a short range much less 
that the interparticle spacing, and, at low temperature, only s-wave scat- 
tering is important. To incorporate these restrictions, it is convenient to 
introduce a simple model, approximating the isotropic interparticle poten- 
tial as U{r) = g5^^\r), where g has the dimension of energy times volume. 
Standard scattering theory for two particles with identical masses M yields 

10 



the identification 

9 = -j^, (16) 

where a is the s-wave scattering length that relates the phase of the scattered 
wave to that of incident wave. For the commonly used bosonic alkali-metal 
atoms (^'^Na, ^^Rb), the scattering length is positive and a few nm (a common 
trapped state of "^Li is unusual in having a large negative scattering length 
that leads to very different physics; see Sec. l4.5p . Typically, the dimensionless 
parameter na^ that characterizes the "diluteness" of the gas is generally 
small. In special cases, however, the scattering length is purposely enhanced 
with an applied magnetic field, as occurs in a "Fano-Feshbach resonance" 
(see Chap. 1). 

3.1. interaction energy for a uniform gas 

In an interacting gas with A^ particles, the ith particle at Vi experiences 
an effective mean-field potential 

Vir,) =gJ25^'\r, - r,) ^ gn{r,), (17) 

where the sum is over all other particles, and the last expression omits small 
corrections of order 1/A^. This effective potential is just the one-body Hartree 
interaction with all the other particles VH(r) = gn{r), familiar from atomic 
and condensed-matter physics. 

It takes an energy E{N+1)—E{N) = gn to add one additional condensate 
particle to a uniform interacting dilute condensed Bose gas in a box. This 
quantity is just the chemical potential /i = {dE /dN)y so that here /i = gn. 
For a uniform dilute gas in a box of volume V , the thermodynamic definition 
of /i implies the total energy EiV, N) = ^gN'^/V. Furthermore, the pressure 
is p = -{dE/dV)N = \gn^. 

3.2. healing length for a uniform gas 

In general, the kinetic energy operator is T = — /i^V^/2M. Since V^ has 
the dimension of an inverse squared length, comparison of the kinetic energy 
T and the Hartree energy Vh = gn leads to a characteristic squared length 

^ 2Mgn 8nna' ^ ' 



11 



where the second form makes use of Eq. (ITB]) . If a uniform Bose gas is 
perturbed locally (by an impurity or a rigid external boundary, for example), 
^ is the length over which the gas heals back to its equilibrium density n. 
In the dilute limit no? <^ 1, the healing length ^ is large compared to the 
interparticle spacing n~^/^^ confirming the mean-field character of the Hartree 
picture of the interaction. 

4. Gross-Pitaevskii picture for a trapped Bose gas 

At zero temperature, a uniform interacting Bose gas in a box has two 
microscopic length scales: the scattering length a and the interparticle spac- 
ing n~^/^ (or equivalently, the healing length ^). In addition, an inter- 
acting trapped Bose gas has a length associated with the oscillator length 
do, which leads to many interesting new effects. Although the original pa- 



pers of Gross [ll| and Pitaevskii [12| focused on a uniform gas, Baym and 



Pethick [13| extended the treatment to a harmonic trap soon after the cre- 



ation of a BEC in 1995 llQl. 

4.1. static behavior 

Gross and Pitaevskii independently started from the time-independent 
Schrodinger equation for the condensate wave function in an ideal gas and 
added the nonlinear Hartree potential, leading to a nonlinear Schrodinger 
equation. With the addition of the trap potential, the Gross-Pitaevskii (GP) 
equation becomes 

^+Kr+^i^r)^=/^^, (19) 

where /i is the chemical potential and \E' is the condensate wave function. 
Strictly, \E' is normalized to the condensate number A'^o, but a dilute Bose 
gas at zero temperature has Nq ^ N so that the condensate normalization 
is generally taken as J dV |\E'p ~ A^. 

For a deeper understanding of this time-independent GP equation, it is 
helpful to start from the corresponding GP energy functional 



EgA% = / dV 



\ "' trap "^ / 

\ kinetic interaction/ 



12 



(20) 



containing the kinetic energy, the trap energy (proportional to the condensate 
density), and the interaction energy (proportional to the condensate density 
squared). Variation of i?Gp[^] — /^-^ with respect to \1/* at fixed normalization 
with the chemical potential as a Lagrange multiplier reproduces the time- 
independent GP equation, Eq. f fT9|) : but note that although this looks like 
a nonlinear eigenvalue problem, the right-hand side of Eq. f lT^ contains the 
Lagrange multiplier fi, rather than the energy per particle. 

For a harmonic potential as in Eq. ( ITTj) . it is convenient to use the mean 
oscillator length do as the length scale and the mean oscillator frequency uo 
times h as the energy scale. The resulting dimensionless form of the energy 
functional in Eq. (120]) has a simple structure, with the first two terms of order 
unity (since they involve only the harmonic oscillator), but the interaction 
term contains a new dimensionless parameter Na/do that characterizes the 
importance of the interaction energy relative to the other two terms. If this 
interaction parameter is small, then the BEC adjusts to balance the kinetic 
energy and the confining trap, like an ideal gas, but the situation is very 
different if this parameter is large. 

Note that the ratio a/do is typically ~ 10~^, but the total number of 
atoms is of order 10^ — 10^, so that this interaction parameter is generally 
large for most experiments. In this case, the repulsive interactions act to 
expand the condensate to a typical mean dimension Rq that is considerably 
larger than d^ (typically Ro/do ~ 10 for Na/do ~ 10'^). Since this expanded 
condensate has a small density gradient, the first term of Eq. (12 Up becomes 
small, and the remaining terms give the approximate Thomas-Fermi (TF) 



energy functional [13 



Etf^ = j dV r^rl^P + \9\'^A , (21) 

which involves only |\l/p and |^l/|^. Variation of this TF energy functional 
with respect to |\l/p immediately yields the TF approximation 

VUr)+g\^{r)\^ = fi, (22) 



which also follows by omitting the kinetic energy in Eq. (1191) . Solution of 
this algebraic equation gives the celebrated TF condensate density profile 

n(r) = ^i^^e[/.-V;,(r)], (23) 

13 



where 6{x) is the unit positive step function. For the general three-dimensional 
harmonic trap in Eq. (1111) . this TF density is an inverted parabola 

where the right-hand side is positive and zero otherwise. Here, n{0) = fi/g 
is the central density, and R'j = 2fj,/Muj^ are the squared condensate radii in 
the three orthogonal directions. 

The normalization condition on the density J dVn{r) ^ N yields A^ 



87m(0)-RQ/15, where Rq = RxRyRz depends on the chemical potential /i 13 
Some analysis yields the important dimensionless relation 

R^n Na 

which is large in the present TF limit. Similarly, the TF chemical potential 
becomes 

/iTF = ^Mu'^Rl = 1^0^, (26) 

showing that /xtf ^ ^o in the TF limit. Note that /^tf is proportional to 
JY2/5. j^gnce the thermodynamic relation ft = OE/dN yields the TF energy 
for the trapped condensate -Etf = f /^tf^^- 

It is natural to use the central density n(0) to define the healing length in 
a nonuniform condensate, and Eq. flTSl) shows that ^^ = dy R^; equivalently, 

do Rq 

The right-hand side is small in the TF limit, so that C, is also small compared 
to the mean oscillator length do. This clear separation of TF length scales 
^ -C (io ^ -Ro is frequently valuable in understanding the physics of the 
trapped TF condensate. 

4-2. time- dependent GP equation 

Equation flT9|) has the intuitive time-dependent generalization (known as 
the time-dependent GP equation) 



'^ dt 



'^ + VUr)+gmr,t)\' 
14 



^(r,t), (28) 



where \l/ now depends on t as well as on r. Comparison with Eq. f ll9p shows 
that a stationary solution has the time dependence exp{—ifj,t/h). To under- 
stand this time-dependence, note that the condensate wave function \l/(r, t) is 
a matrix element of the Heisenberg field operator tpir, t) = e*^*/^^(r, 0)e~*'^*/'*. 
Since the operator ip removes one particle, the states on the right and left 
have N and A^ — 1 particles, respectively. Their energy difference E^ — -Eat-i 
is the chemical potential fi = dE/dN. 

This nonlinear field equation can be recast in an intuitive hydrodynamic 
form by writing the condensate wave function as 

*(r,t) = |^(r,t)|exp[2^(r,t)] (29) 

in terms of the magnitude |\1/| and the phase S. The condensate (particle) 
density is simply n{r,t) = |\E'(r,t)p, and the usual one-body definition of 
the particle current density for the Schrodinger equation shows that j = 
|\l/p^VS'/M. The hydrodynamic definition j = nv then gives the local 
superfluid velocity as 

vir,t) = ^VSir,t) = V^r,t), (30) 

where $ = hS/M is the velocity potential for this irrotational flow. The 
relation between the velocity v and the phase S implies that the circulation 
K = §adl ■ V around any closed path C in the fluid is quantized in units of 



h/M I 



K= — i dl-VS = —ASc, (31) 

M Jc M ' ^ ^ 

because ASc must be an integral multiple of 27r to ensure that the condensate 
wave function \1/ is single valued. 

Substitute Eq. (l29l) into the time-dependent GP equation (!28l) . The imag- 
inary part gives the expected conservation of particles 

oil 

— + V ■ (nv) = 0. (32) 

In contrast the real part gives a generalized Bernoulli equation 
1 h'^ 9$ 



15 



Here, the explicitly quantum-inechanical term involving fi^ is called the quan- 
tum pressure. It is usually small for low-lying collective modes of a trapped 
condensate in the TF limit, since the density is slowly varying. It can be 
important for finite-wavelength disturbances, however, as shown below for 
the Bogoliubov excitation spectrum (see, for example. Sec. 7.2 of 



4-3. Bogoliubov spectrum: linearized hydrodynamic equations for a uniform 
Bose gas 
As a first application of these hydrodynamic equations, consider a uni- 
form stationary Bose gas in a box with periodic boundary conditions (so 
that Vir acts solely to fix the allowed plane waves, but otherwise plays no 
role). Assume small perturbations around the equilibrium number density 
77-0, writing n k. uq + n' , where n' is small. Similarly the velocity potential 
becomes ^q + <^' . The linearized equation of continuity has the simple form 

^ + noV^$' = (34) 

because v' = V<l>'. 

To zero order, the Bernoulli equation yields gno + Md^o/dt = 0. Since 
$0 = hSo/M and gno is simply the equilibrium chemical potential /io for 
a uniform Bose gas, this result reproduces the time- dependent phase of the 
equilibrium condensate wave function \E' oc exp{—ifiot/h), as seen by com- 
paring Eqs (IT9D and (|28|). 

In this example all wavenumbers are relevant, and it is essential to retain 
the quantum pressure in the Bernoulli equation. Thus, it is necessary to write 
^/n ~ y^tTq + ^n'/y/no. The first-order contribution yields the additional 
linear equation 

- ,,, + gn + M—- = 0. 35 

AMno dt ^ ^ 

Assume plane- wave solutions oc expi{k-r—ujt). The two linearized equations 
become 

- iun' - nok"^^' = 0, (36) 

' +g)n' - iuM<l>' = 0. (37) 



\4:Mno 



Solutions exist only if the determinant of coefficients vanishes, which leads 
to the celebrated Bogoliubov spectrum Uk for a dilute Bose gas [ij] (here 
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written in terms of the excitation energy e^ = hwk 



griQh^k^ (hk 



M 



+ 



2M 



2i,2 



gnoh k 

w 



+ i^y 



(38) 



where e° = h'^k'^/2M is the free-particle excitation spectrum. For a summary 
of Bogohubov's original derivation based on quantized field amplitudes, see 
Sec. 7.2 of lij. 

This dispersion relation Uk = ek/h has two distinct regions depending on 
the wave number k. At long wavelengths (small k), the spectrum is linear 
with efc ~ hsk, where s = ygn^jM is the speed of sound. Note that g 
must be positive for this uniform Bose gas, which means that the effective 
interparticle interaction must be repulsive. In contrast, at short wavelengths 
(large fc), the spectrum is quadratic, with e^ ~ e° + gnQ, shifted up from the 
free-particle value by the chemical potential grtQ that reflects the presence of 
the background medium. 




Figure 1: Filled circles show measured excitation spectrum for dilute trapped BEC, with 
solid curve the Bogoliubov spectrum. For comparison dashed curve shows the free-particle 
quadratic spectrum. Reprinted with permission of the authors 15| and the American 
Physical Society. 



.2 



2Mnog/h^ 



The crossover between these two regimes occurs at k^ 
1/^^, where ^ = h/y/2Mnog is the healing length for the uniform gas [see 
Eq. (TT8|) ]. If there were no interactions {g — > 0), then kcr — ^ 0, and the 
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spectrum reduces to that of a free particle for all k. Figure [T] shows the 
measured Bogoliubov excitation spectrum with no free parameters (the free- 
particle spectrum is the dashed line) 

It is interesting to consider the Landau critical velocity ([2]) for this dilute 
interacting Bose gas. Since the Bogoliubov spectrum is linear at small k and 
changes to quadratic at large k, its derivative du^/dk (the group velocity) is 
non-negative. Hence the minimum of e^/hk occurs at long wavelengths (in 
contrast to the case for superfluid ^He) and is simply the speed of sound 



Note that the repulsive interactions {g > 0) are essential to give a nonzero 
critical velocity. For an ideal gas, the spectrum is quadratic for all k and 
the corresponding critical velocity vanishes. Thus an ideal Bose gas at zero 
temperature indeed has a Bose-Einstein condensate, but it is not a superfluid 
in the conventional sense because its critical velocity vanishes. 

4.4- linearized hydrodynamic equations for collective modes of a stationary 
condensate 

A different and very important application of the hydrodynamic equations 
is to study the small low-lying oscillations of an initially static condensate 
with density riQir). In this case, the equation of continuity f l32|) becomes 
dn' /dt + V ■ {riQv') = 0, where n' and v' are small perturbations. Similarly, 
the linearized form of the Bernoulli equation ( 133|) yields gn' + Md^' /dt = 0, 
where V$' = v'. 

The time derivative of the first of these equations and the gradient of the 
second readily yield the generalized wave equation 



which explicitly incorporates the equilibrium condensate density nQ{r) as an 
inhomogeneous effective potential. Note that the quantity gnQ{r)/M can be 
identified as the local squared speed of sound s'^{r). Hence Eq. (HO]) can 
equivalently be written 

^ = V . [s\r)Vn'] , (41) 

which is important in experimental studies of collective modes of trapped 
condensates, as described in Sec. 15.31 



4-5. effect of attractive interactions 

For liquid ^He, the interatomic potential is overall repulsive because of 
the strong repulsive core and the weak attractive van der Waals attraction. 
This is why it remains a liquid down to zero temperature and can only be 
solidified by going to high pressure (many atmospheres) at low temperature. 

The situation is very different for the alkali metals like Li, Na, K, and 
Rb, because each atom is much more polarizable than He. As a result, the 
interatomic potential typically has many bound states, and the sign of the 
s-wave scattering length depends critically on the relative position of the last 
bound state [9|. For example, '''Li atoms have a negative scattering length 
a ~ —1.46 nm equivalent to an attractive interaction. 

As noted in connection with Eq. fl39|) . a uniform dilute Bose gas with 
attractive interactions has an imaginary speed of sound and is thus unstable. 
With a trapped condensate, however, the kinetic trapping energy already 
provides a positive energy that can counteract the negative attractive inter- 
action energy. For a quantitative analysis, consider a spherical condensate 
with oscillator length do and use a trial Gaussian wave function 



rp^ 



where /3 is a variational parameter (note that the condensate radius is Pd^). 
The ground-state energy f l2U]) is easily evaluated to give 



EgiP) = ^Nhwo 



3/1 ^o\ 2 Na 1 



P^ "^ J \ TT do 13^ 



(43) 



and /3 < 1 (/9 > 1) for attractive (repulsive) interactions. Here, the three 
terms represent the kinetic energy, the trap energy, and the interaction en- 
ergy, respectively. This treatment also gives a quantitative justification of 
the approximation made in going from Eq. ( I20i) to ( pTl) . 

It is clear from inspection that Eg{(3) becomes large for /3 — )■ oo because of 
the spatial confinement of the harmonic trap. The detailed behavior for small 
f3 depends crucially on the sign and magnitude of the interaction parameter 
Na/do- For a noninteracting gas with a = 0, the minimum of Eq. ( H3|) 
occurs for /3 = 1, which is just the familiar result for an oscillator. Baym 
and Pethick [l3| used this variational approach to study the expansion of the 
condensate for strong repulsive interactions. The situation is very different 
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for an attractive interaction, no matter how weak, because the (now negative) 
cubic term dominates for /3 — )■ 0, and Eg{(i) diverges to — oo. Hence any 
attractive system is globally unstable with respect to a collapse. For small 
and moderate values of A^|a|/(io, however, Eg{P) does retain a local minimum, 
and the system then remains locally stable. A straightforward analysis shows 
that the minimum disappears at the critical value Nc\a\/do ^ 0.671 and that 
the critical condensate radius is reduced by a factor (3^, = 5~^/^ ~ 0.669. 
For comparison, a numerical study of the GP equation [l6| yields the value 
Nc\a\/do ~ 0.575, which differs from the variational estimate by ~ 17%. 

The attractive interaction in ^Li produces just this behavior. For a trap 
size of order 3 /im, the critical stability limit is Nc ~ 1250 atoms [17]. In 
the experiments, atoms from the surrounding gas continue to augment the 
condensate until it reaches the critical value and then collapses. This process 
repeats in a cyclic fashion. 

5. Selected applications and comparison with experiments 

The GP formalism, either in the form of the time- dependent equation (1251) 
or the pair of coupled hydrodynamic equations fl5^ and fl5^ . is remarkably 
rich and has many applications. Here we consider a few of the most important 
ones. 

5.1. free expansion for different trap aspect ratios 

Much of the experimental data relating to Bose- Einstein condensation in 
trapped gases has been obtained from images of the density distribution of 
an atomic cloud after the trap has been switched off. Indeed such time-of- 
flight expansion of an ultracold gas released from a harmonic trap provided 
some of the first evidence of Bose Einstein condensation [10]. In a series 
of observations taken a fixed time after release the measured size of the ex- 
panded cloud decreases abruptly as the temperature is reduced through the 
phase transition. A second characteristic sign of the transition is observed 
by recording the aspect ratio of the atom cloud (the ratio Rz/R± of axial 
and radial widths) as it expands from an asymmetric trap; for a sufficiently 
long period of free expansion (where the cloud in large compared to its initial 
size) a thermal cloud has an aspect ratio of unity (circular shape) while a 
BEG has an elliptical distribution (as shown in Fig. [2]). In the absence of 
interactions one can explain the difference above as follows. The asymptotic 
density profile of an ideal gas released from a potential refiects the initial 
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Figure 2; A montage of images showing how a Bose-Einstein condensate expands after it 
has been released from a trap and fahs downwards under gravity.The initial elliptical cloud 
is elongated along the horizontal axis and is too small to resolve properly. As it expands 
it becomes circularabout 7 ms) and then elongated vertically. The final shape reflects the 
initial momentum distribution of the BEC wave function and the effects of interactions. 
Data provided by M. Gildemeister and B. Sherlock, Oxford. 



momentum distribution of the cloud in the trap (when the distance pt/M 
travelled by the atoms of momentum p and mass M in time t is much greater 
than the initial size of the cloud) . A cloud of atoms that is far from the quan- 
tum regime has the same velocity distribution (Maxwell-Boltzmann) in all 
directions characterized by the temperature, and hence expands isotropi- 
cally. On the other hand, atoms in a BEC occupy one quantum state and 
are described by a single wave function. The momentum distribution is the 
Fourier transform of the spatial wave function and has an aspect ratio of 
^uJxj'^z when viewed along the ^/-direction. The anisotropy of the confining 
potential is thus reflected in the aspect ratio of the expanding BEC, with the 
aspect ratio being inverted with respect to its initial value in trap, i.e. the 
cloud expands most rapidly in along the axis where the initial confinement 
is tightest, as expected from the Heisenberg uncertainty principle. However 
interactions must be taken into account for a complete description of the 
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expansion of a condensate. They give rise to a mean-field interaction energy 
wliicli is converted to kinetic energy of motion during the release. For larger 
atom numbers this interaction term, proportional to the density, dominates 
the expansion of the condensate (as explained in Sec. 14. ip . Gradients in the 
density drive the initial acceleration of the atoms so the cloud expands more 
rapidly in the most tightly confined direction. This causes the aspect ratio of 
the cloud to invert during expansion as in the non-interacting case, but when 
interactions are dominant the calculation requires solving the set of hydro- 
dynamic equations given in Eq. (132|) an d ([33D , with the external potential Kr 
set to zero at the release time (see also [l3|). This approach based on mean- 
field theory predicted asymptotic aspect ratios for Bose Einstein condensates 
that agree with the experimental results to within a few percent, and this 
was one of the first validations of the mean-field description of interactions 
in dilute atomic vapors. 

5.2. measurement of the condensate fraction 

Superfluid ^He generally forms a macroscopic sample, which facilitates 
the study of the superfluid properties like the normal and superfluid densi- 
ties and persistent currents. In contrast, it has proved extremely difficult to 
measure reliably the low-temperature condensate fraction for bulk uniform 
superfluid ^He. Conceptually, the true many-body ground-state wave func- 
tion \l/o can be expanded in the complete set of states |$j) for an ideal Bose 
gas in a box. The strong repulsive interactions for realistic superfluid ^He 
lead to a linear combination of many ideal-gas states with high energies in 
addition to the noninteracting ground state $o with all the particles in the 
zero-momentum condensate. Correspondingly, the fraction of ^He atoms in 
the zero-momentum condensed state (the overlap |($o|^o)P) is significantly 
reduced, leading to an estimated condensate fraction Nq/N ~ 0.1, in contrast 
to the value 1 for an ideal uniform Bose gas. 

Measurements of Nq/N for superfluid ^He have relied on quasi-elastic 
scattering of high-energy neutrons. In principle, the single-particle momen- 
tum distribution function of bulk superfluid liquid helium has a delta function 
at fc = 0, with a weight proportional to Nq/N. Nevertheless, the experiments 
have proved challenging with considerable uncertainty. Sokol 18| summa- 
rizes the experimental conclusions, confirming the theoretical estimates of 
Nq/N ~ 0.1, and Sec. 8.5 of Ref. [8|] gives an abbreviated account. The 
situation is very different for a dilute ultracold gas. Detection of the con- 
densate fraction is reasonably straightforward and an obvious characteristic 
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feature of the onset BEC in a trap is the sudden appearance of the narrow 
condensate rising out of the broad thermal cloud as T falls below the actual 
transition temperature [ICI]. (Convincing direct measurements of the super- 
fluid properties of quantum degenerate gases followed a few years after the 
first observation of BEC.) 




0.5 1.0 1.5 
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Figure 3: Scaled condensate-fraction ratio No{T)/N compared with that for an ideal 
gas in a general harmonic trap. Solid line is from Eq. ([T5|) with To{N) the transition 
temperature given in Eq. (fTl]). Measured data fit a similar curve but with a reduced 
transition temperature 0.94Tb(A^) (dashed curve). The inset shows the reduction in the 
number of atoms N as the system cools. (Note that N was relatively small in this early 
experiment.) Reprinted with permission of the authors |19| and the American Physical 
Society. 



In the original experiment by the JILA group [19|, and in many subse- 
quent experiments, the atomic cloud was allowed to expand for a definite 
time and then imaged by recording the shadow of the cloud arising from 
absorption of a pulse of resonant laser light. In general, the projected im- 
age has both the broad thermal background and the narrow condensate (for 
T < Tc). The ratio of the areas under these curves yields the condensate 
fraction No{T)/N. Figure [3] shows the measured values which extend down 
to ~ 0.5 Tc. Since the gas is very dilute the condensate fraction closely resem- 
bles that for an ideal trapped gas. Corrections from finite A^ and many-body 
interactions are estimated to be at most a few percent. In recent experi- 
ments 20| the condensation fraction has been measured in a system where 
the scattering length can be tuned to be close to zero so approximating very 
closely to an ideal gas; this improves the ability to observe condensates with 
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small number of atoms and clarifies the interpretation of data like that shown 
in Fig. [31 

5.3. low-lying collective oscillation modes 

Ever since the first observation of BEC in a dilute gas the study of collec- 
tive oscillation modes has been crucial to our understanding of these meso- 
scopic systems. For typical experimental conditions these excitations dom- 
inate the low-frequency response of a weakly-interacting gas to an external 
perturbation. Initial experiments performed at temperatures well below Tc 
tested the mean field theory based on the Gross Pitaevskii equation. The 
very good agreement between theory and experiment with no free fit param- 
eters provided strong evidence for the validity of this model. To model the 
system we return to the linearized hydrodynamic equation Eq. fl4ip . For a 
uniform system the solutions to Eq. (1411) correspond to sound waves prop- 
agating with the speed given in Eq. f l39|) . In a trapped system the wave- 
length of the phonons must be considered in relation to the dimensions of 
the cloud. For wavelengths significantly less than all spatial dimensions of 
the condensate the solutions are sound waves propagating at the local sound 
velocity. Whereas for wavelengths comparable to the size of the condensate, 
the boundary conditions produce standing sound waves at discrete frequen- 
cies which are referred to as the collective oscillation modes of the system. 
Thus the collective oscillation frequencies of a BEC in a harmonic trap do not 
depend on interactions (number of atoms) even though of the speed of sound 
s does, i.e. uj = sk where the wavenumber k is the same for different numbers 
of atoms in the same trap (in the TF regime). This arises because the speed 
of sound is s ~ (riofi'/M)^/^ (see Eq. HO]), and the wavelength of the lowest 
energy excitation is limited by the extent of the cloud in a given direction 
Ri [determined by /x = uqq = MufRf/2 from Eq. f l2B]) in the Thomas- Fermi 
regime]. Hence the lowest mode has a frequency of s/Ri ~ Ui, which is the 
same as for the non-interacting case. 

Collective oscillation modes are classified according to the symmetries of 
the trapping potential. In a spherical potential both angular momentum and 
its axial component are conserved quantities and solutions have the form 
n' = P^ ""^ {r)r^Yim{0, 0) where £ is the orbital angular momentum quantum 
number, m is the magnetic quantum number, P^ "'' (r) are polynomials of 
degree 2nr where rij. is the number of nodes of the radial wave function, 
and Yim{0,4') are spherical harmonics. The corresponding frequencies for a 
spherical harmonic potential are u{nr, i) = uJo{2n'^ + 2nri + Srir + iY^"^ , where 
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loq is the oscillation frequency of atoms in the trap [8|. For rij. = 0, i = 1 the 
motion is a simple dipole oscillation, or sloshing back and forth of the cloud; 
in a harmonic trap this centre-of-mass motion of the many-body system is 
decoupled from its internal degrees of freedom. Thus the dipole modes have 
frequencies exactly equal to those of the oscillations of individual atoms in 
the trap which facilitates calibration of experiments. The modes with orbital 
angular momentum quantum number i = 2 (and n^ = 0) are quadrupole 
modes and in a spherical potential the oscillation frequency is v2ujo for all 
five values of m. The mode with n^ = 1, i = is known as the "breathing 
mode" in which there is expansion and contraction along the radial direction 
changing the volume of the cloud but not its shape and this mode has a 
frequency of a/Swq. Experiments are rarely performed in a spherical trap 
but in many cases the potential does have a cylindrical symmetry about the 
z-axis {ux = Uy = u!±) so that m remains a good quantum number (but not 
i), and two of the three dipole modes remain degenerate. In an axisymmetric 
potential the m = modes with £ = and i = 2 are mixed to give modes 
for which the oscillations along the axial and radial directions are in phase 
(higher frequency) and out of phase (lower frequency), and the \m\ = 2 modes 
remain degenerate (as do the |m| = 1 modes). The frequency spectrum of 
the modes and their shape are shown in Fig. H] for a time averaged orbiting 
potential (TOP) trap with coz/ujr = v8 (see Sec. 12.2 of [SJ for more general 
references) . 

Some of the first experiments on BEG measured the frequencies of the low- 
lying m = mode and an m = 2 quadrupole mode and found good agreement 
with hydrodynamic theory when the number of atoms in the condensate is 
sufficiently large. After this verification, considerable subsequent theoretical 
and experimental effort was concentrated on the so-called scissors modes with 
|m| = 1 as described in Sec. 19.11 The scissors mode differs from the other 
collective oscillation modes because it involves rotational motion about an 
axis (lying in the xy plane) as illustrated in Fig. HI Thus it is sensitive to 
the irrotational flow associated with superfluids. 

Interactions give rise to a nonlinear term in the Schrodinger equation, 
Eq. (IT^ . and this nonlinearity leads to a coupling between the collective 
excitations, e.g. second harmonic generation was observed when the trapping 
potential was adjusted so that the frequency of a high-lying mode was close 



to the second harmonic of the low-lying m = mode [21|. When the BEG 
was driven at the frequency of the lower mode there was a large transfer of 
energy to the mode at twice the frequency. 
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Figure 4: The low energy collective oscillation mode spectrum of a trapped BEC in a 
TOP trap where ojz/ujr = yS'- the three dipole modes are two degenerate modes at oj^ 
and one at Uz', the six quadrupole modes are the ni — low- lying and high- lying modes, 
the doubly degenerate \m\ = 2 modes and the \m\ = 1 or scissors modes. Figure drawn 
by E. Nugent, Oxford. 
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The first experiments on collective oscillations were performed with single 
component Bose gases in three-dimensional harmonic trapping potentials, 
which is by far the simplest situation to analyze. More recently there has 
been work on two-component gases, Fermi gases and in complex and low 
dimensional trapping potentials. The collective oscillation frequencies of a 
BEC do not depend on the strength of the interactions (number of atoms) in 
a harmonic trap for reasons explained above, but they do in a ring trap where 
the lowest-order collective modes are quantized flows around the ring |22| . 
A very special theoretical case is a two-dimensional interacting quantum gas 
confined in a harmonic potential which has a breathing mode at the universal 
frequency of precisely 2uj^, twice the radial oscillation frequency, because of 



underlying symmetry of quantum mechanical Hamiltonian [23[ . A wealth of 
new physics beyond mean field theory occurs in these systems and collective 
oscillations are an important way of uncovering new effects. 

6. Dipolar condensates 

It is helpful to start with a brief review of electric and magnetic dipoles 
and their interactions (Ref. [2J] provides a good general review of these in- 
triguing systems). In SI units, an elementary electric dipole moment p has 
dimensions of charge-length (namely C m) with the electrostatic potential 

where eo is the permittivity of the vacuum. A combination of the dipole 
electrostatic field E{r) = — V$(r) and the energy U = —p ■ E(r) of an 
external dipole in this field yields the interaction energy Udd of two dipoles 
Pi and p2 

J J Pi-p2-3(ri2-pi)(fi2-p2) ..^x 

^- i^^^^TE ' ^''^ 

where ri2 = Ti — r2 is the vector joining the two dipoles. The most notable 
features are the rather complicated angular dependence and the inverse-cube 
dependence on the separation (a long-range interaction, in contrast to the 
usual short-range interaction g8(r) discussed in Sec. 3). 

To simplify this interaction energy, it is common to orient the dipoles 
with an extra applied electric field along 5, in which case Eq. fH5|) has the 
simpler form 

Cdd l-3cos^g 

^'' = 17 rf, ' ^'') 
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where C^d = p'^/^o (assuming identical atomic dipoles) and 9 is the angle 
between the vector ri2 and the z axis. Note the strong angular dependence 
that is proportional to the (negative of the) Legendre polynomial P2{cos6) = 
(3cos^6' — l)/2: two parallel dipoles aligned along the z axis (head-to-tail) 
have negative (attractive) energy, whereas two parallel dipoles side by side 
have a positive (repulsive) energy. It is easy to see that the dimensional 
coupling constant C^d has dimensions of energy- volume, exactly the same as 
the contact interaction constant g from Sec. 3. 

A similar but considerably more intricate argument yields the identical 
expression (146|) for the energy of oriented magnetic dipoles m, where m has 
dimensions of current-area (Am^). The only difference is the definition of the 
magnetic-dipole dimensional coupling constant Cdd = fJ'O^, where /iq is the 
permeability of the vacuum. Note that, by definition, eo/xo = c~^. Compared 
to an elementary electric dipole moment p, an elementary magnetic moment 
m has the same dimension as cp. Furthermore, a typical valence electron in 
an atom has orbital speed of order ac, where a = e^/(47reo/ic) ~ 1/137 is the 
fine-structure constant. The atomic magnetic moment arises from electronic 
motion with velocity v, leading to the typical order-of-magnitude estimate 
m ~ vp, with v/c ~ a. Thus the interaction energy of magnetic dipoles is 



smaller by a factor of a"^ ~ 10"'^ compared to that for electric dipoles 2S ■ 

One approach to parametrize the dipole-dipole interaction is to introduce 
a corresponding dipolar length Odd = CddM/{127ih?). For many purposes, 
the s-wave scattering length a provides the relevant comparison, with the 
dimensionless ratio edd = Odd/a- For ®^Rb, this ratio Odd/a ~ 0.007 is small 
so that dipolar effects are frequently negligible (but see Sec. |8]for an example 
where the dipolar energy is crucial). 

Equivalently, this same dimensionless ratio becomes edd = C'dd/Sf? J24| . 
where g is the contact coupling constant. At present, all experimental studies 
of dipolar effects have relied on the particular atomic species ^^Cr that has a 
large intrinsic magnetic dipole moment m ?» Q^b, where fiB = eh/2me is the 
Bohr magneton. Even in this most favorable case of ^^Cr, the dimensionless 
ratio is only Cdd ~ 0.16. The Stuttgart group has created a BEC of ^^Cr J25|, 
which has allowed many detailed experimental studies. More recently. Lev 
and colleagues have made a BEC of ^^"^Dy, which has an even larger atomic 



magnetic moment /i ~ 10 muB [26 



Recently, powerful laser techriiques have created bosonic heteronuclear 
polar molecules such as ^''K^^Rb 27|, |28| that have large electric dipole mo- 



ments with edd ~ 20 [2J]. The phase-space density required for a quantum 
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degenerate gas of polar molecules has not yet been achieved, but strenuous 
efforts are currently being made in several experiments, and this is an ex- 
tremely promising area for the future work. Such systems will qualitatively 
alter the st udy of ultracold quantum matter and offer many new possibilities. 



Reference [29| provides a recent review of such cold dipolar gases. 

In the context of the Gross-Pitaevskii description, the addition of the 
long-range dipole-dipole interaction potential significantly affects the total 
interaction energy [the third term in Eq. ( 120|) ]. Specifically, the two-particle 
potential now becomes g5{r — r') + Uddij" — f'), leading to the more compli- 
cated contribution to the GP energy functional 

^int[*] = \jdV^*{r)m*{r') [g5{r - r') + U^^{r - r')] ^(r')^(r). (47) 

Here, the first term leads to the local Hartree potential Vh = g\^{r)\'^, but 
the second dipole term is intrinsically nonlocal and couples the total conden- 
sate densities at the two different positions. This second term significantly 
complicates the GP equation, and it is convenient to introduce an additional 
dipole-dipole interaction potential 

Kid(r, t) = j dV U^air - r') \^{r', t)\\ (48) 

that acts in addition to the Hartree potential Vh = S'l^P- The resulting 
time-dependent GP equation 

now contains the nonlocal integral contribution Vdd- For the stationary case, 
the left-hand side of (HOj) simply reduces to /i\E'. 

Based on the form of Eq. (146|) . a cigar-shaped pure dipolar condensate 
should collapse, whereas a disk-shaped pure dipolar condensate should be 
stable. The repulsive contact interaction acts to counteract the collapse, 
and variational studies have confirmed these qualitative ideas, but the de- 



tails become rather complicated [2J]. Assume an axisymmetric trap with a 
particular aspect ratio (this depends on the trap frequencies u^ and u±). Ref- 
erence 30(1 used the BEG of ^^Gr to study the behavior of the condensate for 



various aspect ratios of the external trap. Most importantly, they also used 
a Feshbach resonance to control the s-wave scattering length and confirmed 
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Figure 5: Top row: Experimental images of a dipolar condensate after collapse and ex- 
plosion following a sequence of holding times for release from the trap (the system is 
axisymmetric around the horizontal axis) . Bottom row: Results of numerical simulation 
of the collapse dynamics with no adjustable parameters. The field of view is 130 /im 
X 130 /im. Reprinted with permission of the authors [3l| and the American Physical 
Society. 



the theoretical picture. For example, using a Feshbach resonance to tune 
the scattering length to zero, they demonstrated the stability of a pancake- 
shaped pure dipolar gas. Subsequently, Ref. 3l|] investigated the collapse 
dynamics of the dipolar condensate when the s-wave scattering length is too 
small to support the system. As might be expected from the P2(cos 6) form of 
the oriented dipole-dipole interaction Eq. fH6|) . Fig. [5] shows that the collapse 
dynamics leads to a d-wave symmetric implosion (the system is axisymmetric 
around the horizontal axis). 



7. Mixtures 

This section and the next consider multicomponent BECs, but there are 
two very different types. The first is a mixture of different hyperfine states 
that are confined in a magnetic trap. The second is a "spinor condensate" 
that involves all the magnetic sublevels of a single hyperfine state. Since 
many of the sublevels cannot be magnetically trapped, this system requires 
an optical trap that can confine all the magnetic sublevels. 

For definiteness, we consider the particular species '^^Rb that has been 
used in many of the experiments with mixtures. The nucleus has spin 3/2 
and the single s-state valence electron with spin 1/2 leads to two hyperfine 
states. The lower one has F = 1 and the higher one has F = 2, separated 
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by an energy with an equivalent frequency of 6.8 GHz. 

In a magnetic trap, among the F = 1 manifold, only the state |1, —1) can 
be trapped (using the standard notation \F, mp)) whereas the F = 2 manifold 
has two states |2, 1) and |2,2) that can be trapped. Usually experiments 
focus on the pair |1) = |1,— 1) and |2) = |2, 1) since they have essentially 
the same magnetic moments and thus behave similarly in a magnetic field 
(which means that they experience similar trap potentials). 

7.1. interacting two-component mixtures 

The simplest case of a mixture is two dilute Bose gases with order pa- 
rameters \l/i and ^2 and coupling constants gn and 5^22 for self-interaction 
and gi2 = g2i for mutual interaction. The principal new feature is that 
the GP equation for (say) \&i now includes an additional Hartree term 
VHi2{f) = 5'i2|^2(^)P for the mutual interaction with species 2 
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^i{r,t), 

(50) 
and similarly for the GP equation for \E'2- 

Suppose for simplicity that the two condensates are uniform. It is not 
difficult to see that the conditions for miscibility are [9|, Sec. 12.1] 

9n > 0, 5-22 > 0, and 5-115-22 > 9i2- (51) 

The first and second conditions ensure that each pure condensate is sta- 
ble against collapse (as is familiar for a single- component condensate). In 
contrast, the third condition ensures that the overlapping mixture is stable 
against phase separation. If this latter condition is violated, then the uniform 
two-component system will phase separate, and the long-wavelength collec- 
tive modes of the coupled uniform system will have imaginary frequencies. 

The situation in a trap is somewhat different, since the harmonic trap 
acts to confine the two components. Unless the self-interaction constants are 
essentially equal, one component will form a denser core, surrounded by (and 
partially overlapping with) the other less-dense component that experiences 
the inward pressure of the trap and outward pressure of the inner core (anal- 
ogous to density stratification in an inhomogeneous self-gravitating body). 

As explained below, external electromagnetic fields allow a convenient 
transformation between these hyperfine species, for example a sudden tran- 
sition from |1) to |2). Reference [32!] modeled this behavior with a single- 
component time- dependent GP equation with a time- dependent scattering 
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length a(t) that changed from an to 022 discontinuously at t = 0. After 
a variable time-delay, the trap was turned off, allowing the condensate to 
expand. Initially, the condensate shrank, followed by compressional oscilla- 
tions. This behavior is understandable because an > 022 for these particular 
states. Thus the initial density profile for |1) could no longer sustain the 
equilibrium radius once the species transformed to |2), experiencing radial 
shrinkage and radial oscillations. A fit to the observations indicated that 
011/022 ~ 1-06. 

7.2. electromagnetic coupling between two hyperfine states 

The previous model in fl50p assumes that the two components inter- 
act only through the density (the mean-field Hartree terms) that contain 
no phase information. Specifically, each condensate has the representation 
^i(^) ~ l^i(^)l 6^P[^'S'j(r)]. Each order parameter is single- valued when r 
executes a closed path at any given instant, so that the phase is 27r-periodic 
for any such closed path. As noted in connection with superfluid ^He, the cir- 
culation for each component is therefore quantized in units of h/M. Techni- 
cally, each complex one-component order parameter has U{1) symmetry. Far 
below the superfiuid transition temperature, the quantized circulation "ex- 
plains" the existence of persistent currents because any fiuctuation that can 
change the circulation is very improbable, reflecting the topological charge 
associated with the 2it periodicity. 

The situation is very different when external electromagnetic flelds cou- 
ple the two components, leading to off-diagonal terms in the combined GP 
equations. The resulting total order parameter no longer has two separate 
complex scalar functions [each with f/(l) symmetry]; instead it now has two 
components with intrinsic coupling, qualitatively changing its character into 
a single SU{2) structure. For elementary discussions of this essential model 



system, see, for example, J33|, Chaps. 7-11] and j34. Sec. 4.4]. The dynamics 
of this spin- 1 becomes similar to that familiar from nuclear magnetic reso- 
nance, obeying what are here called the "optical Bloch equations." In the 
present case, the transformation from |1) to |2) is analogous to a vr pulse 
that rotates a spin up into a spin down. If the pulse is twice as long, the 
resulting 27r rotation reproduces the initial |1) state apart from an overall 



phase. Experiments [32| verify the periodic transfer between the two states 
as the length of the apphed pulse varies (each hyperfine state has a shghtly 
distinct resonant frequency, so that each can be imaged separately). 



32 



This electromagnetic coupling has the remarkable feature of allowing a 
direct measurement of the relative phase between the two components, as 



becomes clear from the dynamical equations of the coupled system [35 



n— (^^\ ~ /^'^ + ^*^i + ^^1 + ^^12 \ hn{t) exp(2a;rft) \ f^i 
dt \^2) ^ V I ^^(^) exp(-^u;rft) T + V,,2 + Vh2 + VH12) \^2 

(52) 
where T denotes the kinetic energy, u^i is the frequency separating the two 
hyperfine states, and VL{t) reflects the two-photon electromagnetic coupling 
between the two hyperfine states. It is known as the "Rabi" frequency and 
depends on the strength of the electromagnetic coupling field (the time de- 
pendence here arises because it can be turned on and off, but more general 
situations can and do arise). Transforming to the rotating frame eliminates 
the explicit rf time dependence. 

The experiment starts with population in |1) and applies a 7r/2 pulse that 
creates a linear superposition of the two states, then waits a variable time and 
applies a second 7r/2 pulse. The trap is turned off, expanding the condensate. 
Resulting images of either component show oscillations depending on the 
length of the variable delay, which provide a measure of the relative phase. 

This and other experiments indicate that the electromagnetic coupling 
has the remarkable property of changing the topology of this two-component 
system. When the coupling is off, they are simply separate U{1) complex 
scalar order parameters, each with its own phase angle and corresponding 
quantized circulation (each has the topology of a cylinder). When the cou- 
pling is on, they instead form a single coupled SU{2) system with two angles 
on the polar sphere and no topological quantization (the topology is now that 
of a sphere). The JILA group studied this system in considerable detail, and 
its topology proved essential in the first experimental creation of a quantized 
vortex in a dilute BEC [36], as discussed in Sec. El 

8. Spinor condensates 

The previous section focused on mixtures that usually involve two dif- 
ferent F hyperfine manifolds in a magnetic trap. In this situation, only the 
weak- field seeking states are relevant, since the other states no longer remain 
confined. In contrast, a spinor condensate involves a single hyperfine mani- 
fold, but it uses an optical trap that retains all the mp magnetic states (the 
simplest and usual case is F = 1 so that the mp levels are +1, 0, and —1). 
See Chap. 1 for an introductory discussion of these traps. 
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8.1. spinor condensates: special case of F = 1 

The possibility of trapping all three magnetic sublevels in an optical trap 
rapidly led to the study of spin-1 Bose-Einstein condensates (as well as more 



general cases) [38|,l39[. In contrast to the previous mixture of distinct species, 
the rotational invariance of the interactions between two spin-1 atoms now 
leads to special restrictions. Here, the macroscopic order parameter is a 
three-component vector (written as a transpose for convenience) 

*^ = (^1 ^0 ^~i) (53) 

involving the three mp states. Only s-wave scattering is relevant in the 
present low-energy limit, and the interaction between the atoms has the 
familiar contact form 

VUri - r,) = 6{r, - r,) J] ^^ P^, (54) 

F 

where F is the magnitude of the total hyperfine spin of the two atoms F = 
Fi + F2 and Vp is the projection operator onto the appropriate value of the 
total F. 

For bosons each with hyperfine spin Fj = 1, the only allowed total values 
are F = and F = 2, with two scattering lengths Qq and a2- The interaction 
potential can then be rewritten in the equivalent form 

V^int(ri - rs) = 5{n - ^2) (^0 + 92 F^ ■ F^) , (55) 

where the effective interaction constants are given by 

47r^^ 2a2 + ao A-nh^ 02 - ao ,^^. 

In practice, the two scattering lengths are roughly comparable, so that \g2\ <^ 
5^0- Specifically, g2 is small and positive for ^^Na, but it is small and negative 
for ^^Rb, which has crucial practical consequences for the distinct behavior 
of these two atomic species. 



Ho [38[ introduced an effective energy functional, writing the spinor order 
parameter in the form '^ai'f) = \/n{r) Ca{f), where the spinor index runs 
over the values 1, 0, —1, n{r) is the common density for all three components, 
and C«(^) is a normalized spinor with (^'1' ■ ^ = 1. The ground state follows 
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by minimizing the energy with fixed total particle number (enforced with a 
chemical potential), leading to an effective energy functional K = E — fiN 



K = dV 



^' '-H' + ^\^Cf -n[f^-VUr)] + ^{9o + g2{Ff) 



2M I ' I 2M ' '' '^ -V /J ' 2 

(57) 
Here the trap is taken as independent of the hyperfine state, and fi determines 
the total number of atoms. 

Apart from the gradient terms, the ground-state spinor („ follows by 
minimizing the spin-dependent part of the energy ^n'^g2{F)^, where (F) = 
Sofl CQ-^a/3C/9 is the expectation value of the appropriate spin-1 matrices |9|, 
Sec. 12.2.1]. If g2 is positive (as for ^^Na), then the minimum occurs for 
I (F) I = 0. Such states are called "polar" and they are those obtained by spa- 
tial rotations of the hyperfine state \mp = 0). If g2 is negative (as for '^^Rb), 
then the minimum occurs by maximizing the average spin with \{F)\ = 1. 
These states are called "ferromagnetic" and they are those obtained by spa- 
tial rotation of the hyperfine state Imp = 1). 

8.2. experimental studies of spinor condensates 

Soon after the optical trapping of a BEC of ^^Na atoms, the MIT group 
made a detailed experimental study of these fascinating spin-1 spinor sys- 
tems, where the spin-dependent interaction constant g2 in (l56ll is small and 
positive, favoring polar configurations with (F) = 0. For definiteness, we 
focus on the formation of ground-state spin domains in an external magnetic 
field [40]. The resulting spin structures can be either miscible or immiscible 
depending on the applied magnetic field and on which of the nip components 
are occupied. They imaged the spin domains by time-of-fiight expansion fol- 
lowed by a Stern-Gerlach separation of the various spin components in an 



inhomogeneous magnetic field. Reference J4l| comprehensively reviews this 
and related experiments on the polar spinor condensate ^^Na. 

Recently, Stamper-Kurn's group at Berkeley has studied the very different 
situation in ®^Rb, where the relevant spin-dependent coupling constant (?2 
in (13^ is small and negative, favoring ferromagnetic spinor structures with 



\{F)\ = 1 |42|. Here, the relevant scattering length is not the usual s-wave 
a, but the much smaller spin-dependent quantity Aa = (02 — ao)/3 that 
follows from Eq. fl56|) . Consequently, the usually small dipolar length a^d 



now becomes crucial. Specifically, the dimensionless ratio |add/Aa| ^ 0.4 is 
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no longer small, which leads to many intriguing new phenomena involving 
the dipolar energy (effectively a dipolar quantum ffuid |42|). 
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Figure 6: Spontaneous dissolution of helical textures in a quantum degenerate ^^Rb spinor 
Bose gas. A magnetic field gradient prepares transversely magnetized (b) uniform or (a), 
(c) helical magnetic textures. The transverse magnetization after a variable time of free 
evolution is imaged in the xz plane, with orientation indicated by the hue and amplitude 
brightness, (b) Uniform texture remains homogeneous for long evolution times, whereas 
(c) with initial pitch A = 60 /zm dissolves over ~200 ms, leading to a sharply modulated 
texture. Reprinted with permission of the authors 42] and the American Physical Society. 



The experiment involved an asymmetric triaxial trap with a thin elon- 
gated elliptical two-dimensional TF condensate that initially had a helical 
spin texture with a variable pitch (periodicity length A between 50 and 150 
/im). This helical spin texture evolved for a variable time, and the vector 
magnetization was then measured nondestructively. For comparison, they 
also prepared similar uniform samples with much larger pitch (A S> Rz)- 
Figure [6] shows the evolution of these initial textures (a), including both (b) 
uniform and (c) helical structures. For (c), note the small-scale structure 
of the magnetization of order 10 /im, much smaller than the initial helical 
pitch. Several experimental tests associate this short-range modulation with 
the dipolar interaction. 
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9. Rotating Bose gases and quantized vortices 

The superfluid properties of liquid lielium at low temperatures (below 
the lambda-point) were apparent soon after liquefaction was achieved at the 
beginning of the 20th Century, in particular the absence of viscosity that 
allows the liquid to leak through small channels. The properties of superfluid 
helium were explained in an intuitive phenomenological way by the two- 
fluid model (of Landau and others) outhned in Sec. II. Ij although there is 
a connection with Bose-Einstein condensation, it has not been central to 
understanding the superfluid helium as can be seen by comparing the physics 
discussed in Sec. I1.3l with that in 14.31 

The study of BEC in dilute atomic gases developed in a very different 
way. The first experimental breakthroughs confirmed the prediction of Bose 
and Einstein and found that the critical temperature was that predicted by 
statistical mechanics (within a few percent). About four years after the first 
BEC was created, experiments started to demonstrate the striking superfluid 
properties that distinguish it from a classical fluid. This section mainly 
concerns those methods that use rotation to probe the quantum properties. 
For completeness, however, we mention the work on the flow of the quantum 
fluid around a strongly repulsive potential created by a focused blue-detuned 
laser beam. Displacing this laser beam back and forth acts like a 'macroscopic 
object' moving inside the cloud of atoms; it was found that there is a critical 
velocity [defined in Eq. ( 139|) ] below which there is no observable dissipation 
of energy. 

9.1. the scissors mode 

Some of the first evidence for the superfluidity of dilute atomic gases in 
the BEC regime came from the quantitative measurements of the so-called 
scissors mode which derives its name from nuclear physics where there is 
an excitation of the protons and neutrons in deformed nuclei (as explained 
below) |43[. In the case of BEC the cloud of ultracold atoms is non-spherical 
because of asymmetry in the trapping potential, i.e. the three frequencies in 
Eq. (ITT]) are not all the same (which is almost always true in practice), so 
that the density contours are deformed into ellipses. When the potential is 
suddenly rotated through a small angle about a symmetry axis, the deformed 



density distribution finds itself offset from its equilibrium position [4^ and 
therefore starts rotating to and fro through a small angle without changing 
shape corresponding to the modes labelled \m = 1| in Fig. |H Superfluids 
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behave in a distinctive way in rotating systems because of their irrotation 
flow (V X fs = as discussed in Sec. Il.ip : this leads to a reduced moment 
of inertia as compared to the classical value for a 'solid-body' with the same 
density distribution. 

In trapped three-dimensional ultracold gases the superfluid and BEC frac- 
tion are considered to be the same (the condensate fraction is small in helium, 
as discussed in Sec. 15. 2p . Two-dimensional systems, however, can be com- 
pletely superfluid even when there is no BEC (see Sec. l9.2p . and there is great 
interest in investigating the relationship between superfluidity and BEC. 

As mentioned above the scissors mode experiments with ultracold gases 



were suggested j43[ based on an analogy with groundbreaking work on su- 
perfluidity of nuclear matter. The scissors mode in strongly deformed nuclei 
was first suggested in a model where the protons and neutrons are assumed 
to be two distinct deformed quantum fluids, e.g. two ellipsoidal distributions 
with a common center. These deformed distributions rotate relative to each 
other around a common axis (two-rotor model) like the blades of a pair of 
scissors about the pivot. 

9.2. the nucleation of vortices 

Perhaps the most obvious way to create a vortex in a fluid is rotating 
the container as in experiments with superfluid helium, or using laser beams 
moving with a circular motion to stir the system like a spoon in a cup of tea. 
However the first vortex created in a BEC of ultracold Rb atoms was achieved 
at JILA by an ingenious scheme which directly controlled the phase of wave 
function, in a manner far removed from any previous work (on helium) J45| . 
That experiment used a mixture of ultracold atoms in two different magnetic 
states as described in Sec. [71 By driving a two-photon transition between 
the states (with a combination of microwave and RF radiation), a ring of 
atoms was created which had a phase winding of 27r corresponding to a vor- 
tex state, like a persistent flow with one unit of angular momentum (and 
further work on persistent flow is discussed below). In this scheme the cen- 
ter of the ring was filled with atoms in a different Zeeman sub-state, but 
these atoms were selectively removed to leave a cloud of atoms containing a 
more conventional vortex (and much of the rich physics of vortices and other 
topological structures in multi-component systems remains to be explored). 

In principle this scheme is very flexible, but much of the subsequent work 
has used experimental methods in which mechanical rotation or stirring im- 
parts angular moment to the condensate. Analogous to the "rotating bucket" 
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experiment in superfluid ^He, ultracold atoms have been trapped in a rotat- 
ing potential whose contours of constant energy are elhptical (at a given in- 
stant of time). However once a BEC has formed and estabhshed long-range 
phase coherence there is significant hysteresis: the condensate can remain 
in a (metastable) state without vortices even when the confining potential 
rotates sufficiently fast that vortex states have lower energy. The nucle- 
ation of vortices requires a way to overcome the energy barrier. This has 
been achieved by driving the system at a frequency that resonantly excites a 
quadrupole mode so that deformation builds up until the cloud of atoms is 
sufficiently perturbed that vortices can enter; vortices form in the low-density 
outer reaches of the cloud and then move towards the center. 

The characteristic core size in the healing length ^ given in Eq. fllSp . 
The wave function for a vortex has the form \l/(r, 0) = f{r/^)e^'^, and a 
variational solution of the GP equation gives the approximate amplitude as 
f{x) = x/y/x^ + 2 which provides a good fit at both large and small x = r/^ 
(when compared to more accurate numerical solutions). 




Figure 7: Images of expanded Bose-Einstein condensates of ®^Rb atoms showing (a) small 
vortex array for slow rotation, and (b) large vortex array for rapid rotation. Note the 
highly regular triangular form, (c) An array of vortices that has not settled down to a 
regular lattice as in the previous images; stirring or shaking such a system creates a tangle 
of vortex lines corresponding to a turbulent quantum fluid. Images (a) and (b) reprinted 
with permission of the authors |46l| an d the American Physical Society. Image (c) provided 
by R. WiUiams and S. Al- Assam [47} . 

Rather than starting with a static BEC and trying to break up the es- 
tablished phase, an alternative strategy for nucleating vortices is to spin up 
a thermal cloud of trapped atoms and then cool them below the critical tem- 
perature in the rotating frame. Arranging preferentially to remove atoms 
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with less than average angular momentum during evaporation increases the 
rotation rate as the temperature decreases. This puts the system directly 
into the equilibrium state for the given rotation rate, i.e. the quantum fluid 
has the number of vortices corresponding to the rotation rate of an equivalent 
classical fluid, see Fig. [7] (as outlined in Sec. II. 2p . 

Generally the vortices are too small to be observed in a trapped BEG and 
images such as Fig. [7] are taken after a time-of-flight expansion. Indeed for 
pancake-shaped condensates the vortices expand proportionally more than 
the overall size of the cloud to give extremely high contrast images. By 
observing along the axis of rotation the cores of the vortices appear as dips 
in the density of the cloud of atoms (in images recorded by absorption of 
resonant laser light). 

Recently another method has been developed that is closely related to 



the original phase imprinting scheme j47j. A Bose- Einstein condensate is 
loaded into a deep 2-D optical lattice to give an array small condensates 
(containing hundreds of atoms at each site); the phases of these separated 
groups of atoms can evolve independently so that when the optical lattice 
potential is rotated, about an axis perpendicular to the plane of the lattice, 
phase differences build up corresponding to the velocity field in a rotating 
system [see Eq. (!30|) ]. Thus when the lattice depth is reduced adiabatically 
to zero the small clouds join together with relative phases that give vortices 
throughout the system (with no hysteresis). Vortices rotating in the same 
direction repel each other and they arrange themselves into a lowest energy 
configuration that is a regular triangular array (or Abrikosov lattice similar 
to that of vortices in superconductors), assuming that there is some damping 
in the system. 

In some situations, the regular arrangement of vortices aligned along the 
rotation axis as in Figs. [7] can be disrupted to create a turbulent system. In 
superfiuid ^He, the much smaller healing length and larger sample size allows 
for many vortices and quantum turbulence in the evolution of a tangle of 
vortices has been studied. Another fascinating type of system is ring vortices 
that resemble smoke rings in the condensate; these have been created in 
experiments with ultracold atoms by producing a sharp phase discontinuity 
in the BEG that evolves into a ring vortex (albeit with less contrast than 
for vortices whose cores go right though the cloud). In addition to studying 
the properties of individual vortices, these systems allow the study of the 
collective properties of "vortex matter". Interesting effects are predicted to 
occur when the number of vortices becomes very large, corresponding to very 
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fast rotation, as described in the following section. 

In addition to the work on BEC, vortices have been used to investigate 



superfluidity in quantum degenerate Fermi gases [48|, which has connec- 
tions with other systems such as superconductors, neutron stars etc. Cold 
atomic vapors have the important characteristic that the strength of inter- 
actions between atoms can be precisely controlled over an enormous range 
(for both fermions, bosons and mixtures) thus making these systems a fan- 
tastic test bed for theoretical ideas. The experiment of Zwierlein et al. at 
MIT 48| used the fermionic isotope of lithium (^Li). Adjusting the mag- 



netic field experienced by the atoms to increase the strength of interactions 
(near a Fano-Feshbach resonance where the scattering length is resonantly en- 
hanced) causes pairs of atoms to bind together to form long-range molecules, 
i.e. molecules in a very high-lying vibrational level of the molecular poten- 
tial where the internuclear spacing is orders of magnitude greater than that 
of a Li2 molecule in its ground vibrational level. These molecules are com- 
posite bosons and hence at low temperature they are in a BEC. A vortex 
lattice was created in this BEC by imparting angular momentum to the 
system by rotation, as in experiments with bosonic atoms. The binding en- 
ergy was decreased (by changing the magnetic field) so that the ultracold 
molecules became weakly bound pairs of fermions closely analogous to the 
Cooper pairs (of electrons) in superconductors. This allowed a detailed study 
of the crossover from a BEC of molecules to a Bardeen-Cooper-Schrieffer su- 
perfiuid of loosely bound pairs, and this particular experiment showed that 
superfluidity as evidenced by the existence of vortices persisted in this BEC- 
BCS cross-over. The ability to scan the system from one regime to another 
is a truly remarkable feature of ultracold vapors. 

The first observation of a vortex in an ultracold gas was a challenging 
experiment but nowadays vortices have been observed in a great variety of 
ways; even in non-rotating systems where vortices can (sometimes) arise 
spontaneously after a rapid quench, i.e. sudden cooling through the phase 
transition, or by merging Bose-Einstein condensed clouds that have different 
(random) phases. Another method that has be successfully used to nucleate 
vortices spontaneously, without rotation or stirring, is by merging of multiple 
trapped BECs |49|. In this experiment three BECs were formed that were 
separated from one another and then the potential barriers between them 
were lowered so that they merged together. The independent, uncorrelated 
BECs had different (random) phases and sometimes there was sufficient phase 
winding for vortices to form. 
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The physics of the persistent flow of a superfluid around a ring and a 
quantum vortex are closely linked. For a ring-shaped cloud of BEC the 
atoms are at a fixed radius and therefore the condition on circulation in 
Eq. fl3T|) corresponds exactly to quantisation of angular moment in units of 
h. Note, however, that vortices with more than one unit of circulation are 
energetically unstable, e.g. a vortex with two units rapidly breaks up into two 
vortices. This has been used to measure angular momentum of flow around a 
ring formed by the combination of a repulsive barrier (from a laser beam with 
blue frequency detuning) passing through the center of a trapping potential; 
slowly reducing the height of the central barrier to zero puts the atoms into a 
harmonic trap where the ring-shaped cloud rapidly separates in many vortices 
of unit circulation which can be counted. Vortices with circulation greater 
than unity can be energetically stable in anharmonic potentials, e.g. when 
there is a positive quartic (r^) term in addition to harmonic confinement (r^). 
Such traps are of interest in the pursuit of the quantum states that arise at 
very fast rotation rates as described below. 



9. 3. the use of rotation for direct quantum simulation 

The rotation of neutral atoms can be considered as a way of simulating 
the effect of a magnetic field B on charged particles. A simple illustration 
of this analogy can be seen by comparing the form of the Lorentz force on 
a charged particle in a magnetic field, q{v x B), and the Coriolis force on a 
neutral particle of mass M in the rotating frame, 2M{v x $7), suggesting a 
mapping 

qB = 2Mn. (58) 

where we generally take 17 = Qz. More rigorously one can consider the form 
of the Hamiltonian of a single neutral particle in a 3D harmonic trapping 
potential in the rotating frame, 

Hn = Ho-n-L (59) 

= ^[P- M{n X r)f + lM(a;i - ^^){x^ + V^) + \Mulz\m 

It is apparent that the kinetic energy term is formally equivalent to the 
usual gauge-invariant term (jp — A)^/2M, with S = V x A, that arises for 
a charged particle in a magnetic field. (We take the effective charge as 1 
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so A has the dimension of momentum [50[.) Note that the centripetal ac- 
celeration which arises from rotation reduces the radial harmonic trapping 
frequency, uj± — >■ (cu^ — ^2^)^/^. Thus in experimental conditions where the ro- 
tation frequency fi tends to the radial trapping frequency cj± radial trapping 
becomes very weak and the atomic cloud expands to be a very thin pancake- 
shape. Under these conditions the temperature of the quantum degenerate 
gas is so low that ksT <^ huJz and the motion in this direction is frozen 
out. This makes the Hamiltonian of Eq. (!60|) almost equivalent to that of 
a two-dimensional electron gas under a magnetic field (but with some resid- 
ual harmonic potential in the radial direction). Interactions have not been 
included in this description but typically the chemical potential /i = gn{0) 
is small in this regime because of the cloud has low density, n(0), when it 
has spread out radially. The close analogy between the physics of rotating 
neutral atoms and electrons under a magnetic field has led to considerable in- 
terest in the possibility of achieving strongly-correlated quantum Hall states 



in a rapidly rotating atomic gas (for a recent review see |51|). 

The use of rotation for simulating a gauge field with ultracold atoms has a 
significant drawback, namely exotic strongly-correlated states are predicted 
to occur in a Bose gas rotating very close to the centrifugal limit, Q ~ u±, 
at which rotating atoms are expelled from the trap. There are considerable 
technical challenges associated with rotating close to the centrifugal limit, e.g. 
the harmonic trap must be very axisymmetric {ux very close to Uy), and any 
static perturbations in the lab frame lead to substantial heating. These issues 
have prevented fractional quantum Hall effect physics from being achieved 
in a rotating Bose gas, though impressive attempts have been made [52|. In 
order to circumvent the pragmatic problems associated with rotation there 
has been a development of alternative methods for creating synthetic gauge 



fields for ultracold atoms [50[, as described in the following section. 



10. Synthetic (artificial) gauge fields — vortices without rotation 

This particular vector potential A = Mflxr in Eq. (!60|) yields an effective 
magnetic field Sgg = V x A = 2Mf2 that is uniform and the corresponding 
vector potential is analogous to what is called symmetric gauge. (The crucial 
requirement for a nonzero magnetic field is that the vector direction of A 
differs from that of the intrinsic spatial dependence: for symmetric gauge, A 
is along (p and its magnitude depends on the distance r.) Other choices yield 
the same magnetic field if they are related through a gauge transformation 
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A — ^ A' = A + VA, where A is a scalar gauge function. Under such 
a gauge transformation, the quantum-mechanical state vector undergoes a 
phase change -ip ^ ip' = exp(iA/h)^ J3J, p. 200]. For example, another 
common choice takes A' = —2MQ y x (Landau gauge) where the constant 
coefficient is chosen so that the two vector potentials yield the same effective 
uniform magnetic field V x A' = 2Mf2. 

The line integral of A plays an essential role; the special choice of a closed 
contour C gives the familiar gauge- invariant result §^A-dl = J dS ■ V x A, 
which is just the magnetic fiux enclosed by the contour. In the present case of 
unit fictitious charge, the fiux quantum is h, and the enclosed fiux is simply 
h X A^„, where N^ is the number of fiux quanta enclosed by C. Equivalently, 
the enclosed fiux is h times the net phase change around the closed contour 
(namely h x 27rN^). In particular, if a quantum-mechanical state yields a net 
phase change when it executes a closed path C, it experiences a gauge field 
A (either from real magnetism or from an artificial/synthetic gauge field). 

This observation leads to the important concept of geometric phase, often 



called Berry's phase [53|, |3J, Sec. 10.2]. If a nondegenerate quantum eigen- 
state |xi(^)) depends on r along a contour C, then there is a net geometric 
phase change on once encircling the contour C 



'yc = h-' idl-A{r), (61) 

where the integrand 

Air)=th{xiir)\Vxiir)) (62) 

is a synthetic or artificial vector potential that acts just like any familiar 
applied gauge field [50|, Sec. LA]. The fine integral 7c is real and can be 
interpreted as the Berry's phase |53J acquired by an atom after adiabatic 
transport around the closed loop C in coordinate space. This quantity is 
called geometric phase because it depends only on the path and not on the 
speed (assumed slow). 

To make a connection with cold atoms, it is helpful to summarize a simple 
toy model used in [501, Sec. I]. Specifically, consider a two-state atom with 
bare states \g) and |e), and a one-body Hamiltonian p^/2M + Vtr that is 
diagonal in this basis. Add a couphng operator U that refiects, for example, 
electromagnetic coupling between the two bare states. Assume the general 
hermitian traceless form 

2 Ve^-^sin^ -cos^ ) ' ^ > 

AA 



where ^2 is a generalized Rabi frequency [compare Eq. fl52p ]. is a mixing 
angle and determines the phase of the off-diagonal coupling terms. The 
2x2 matrix is the product n ■ cr where the unit vector n is characterized 
by the spherical polar angles 6, 0, and cr = ^,=i Xjaj represents the Pauli 
matrices. All three parameters Q, 6, and (p can depend on position r. 
The operator U has two eigenvectors 

I . / cos(^/2) \ I . /-e-*'^sin(^/2)\ ._., 

I^^^=^e^^sin(^/2)J 1^^^=^ cos(^/2) ) ^^^^ 

with eigenvalues hQ/2 and —Ml/2, respectively; these eigenvectors are usu- 
ally known as dressed states and are the eigenstates of the atoms in the 
radiation field. They can serve as the basis for an expansion of the full state 
vector of the particle: 

2 

|M/(r,t)) = 5^^,(r,t)|x,(r)). (65) 

i=i 

Consider the action of the momentum operator p = —iK^ on the state 
\^{r,t)). Since the state vector fl65|) is a product of two position-dependent 
functions, this operation produces p\^) = J2j=i[{P'^j)\Xj) + i^jPlXj)]- A 
little manipulation yields 

2 

p\^) = Y^ (P^Jk - ^jk) ^k\Xj) (66) 

j,k=i 

with Ajfc = ih{xj\'^Xk) a 2 X 2 vector matrix in the dressed-state indices, 
and we have inserted the 2x2 completeness relation to obtain the second 
term in Eq. f lBB]) . 

If the system is prepared with the particle in state |xi), and the particle 
moves sufficiently slowly, it will adiabatically remain in this state. The pro- 
jection of the kinetic-energy term (p^/2M)|\E') onto the state |xi) yields the 
effective Schrodinger equation for the amplitude ipi describing the center-of- 
mass motion in the first internal state 



dipi 



^^-^^' 'V. + '^ + W 



2M 2 



^1, (67) 



where the induced vector potential A is given in Eq. ( l62|l and W = \Ai2\^ /2M 
is an induced scalar potential arising from the elimination of the second 
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atomic state \x2)- Note that this synthetic gauge field A means that cold 
neutral atoms can serve to study magnetic effects usually associated with 
electromagnetic charge. 

In the present case, the synthetic vector potential has the form A{r) = 
^h{cos9 — 1)V0, and the effective magnetic field becomes B{r) = V x A = 
I^V(cos^) X V0. For a nonzero effective magnetic field, V^ and V(/) must 
not be collinear. 

Among the many proposals for experimental implementation of these 
ideas 50|, Sec. II], we here focus on one scheme 5J| that has proved prac- 



tical in creating vortices with a synthetic gauge potential and no externally 
applied rotation |55|. Two slightly detuned laser beams counter-propagate 
along X with a difference wave vector qx. Together these beams induce 
Raman transitions, and the corresponding Rabi frequency has a spatial de- 
pendence oc e'^'^^ that appears as the phase (p = qx oi the coupling matrix 
U in Eq. f l63|) . In addition, a uniform magnetic Zeeman field along y pro- 
vides detuning away from the resonant coupling and determines the mixing 



angle 9. Detailed analysis in momentum space for this two-level system [54 
shows that the minimum in the dispersion relation shifts away from k^ = 
to a finite value with an approximate dispersion relation {hkx — AxY /2M. 
(Reference |50|, Sec. II. D] gives the corresponding real-space analysis.) This 
shift yields what is effectively a synthetic vector potential A^ that depends 
on q and on the Zeeman field. A^ would be uniform in the absence of ad- 
ditional perturbations. In the presence of a time-dependent Ax{t), however 



the derivative —dAx/dt induces an effective electric field E^. Reference [56 
ramps a uniform A from an initial value to a final value and measures the 
resulting pulsed electric field, obtaining values for both the canonical mo- 
mentum p and the mechanical momentum Mv = p — A. 

Application of a variable magnetic field gradient along y supplies the nec- 
essary spatial variation Ax{y) like the Landau gauge, and the combination 
of 0(x) and 6{y) together produce an effectively uniform rotation (synthetic 
magnetic field) over a reasonably broad spatial range. In practice the ac- 
tual experiment involves all three Imp) = 1, 0, —1 magnetic substates, and 



the theory requires numerical analysis. Nevertheless, Ref. [55[ was able to 
produce vortices without rotation, as shown in Fig. [HI One puzzle is that 
the vortices do not form a regular array, in contrast to general theoretical 
predictions, but this may arise from insufficient time to reach equilibrium. 

This field is developing rapidly, especially in connection with trapped 
atoms in optical lattices (see Chaps. 1 and 5 for general discussions). Refer- 
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Figure 8: Vortex formation for increasing applied magnetic field gradient with N — 1.4 x 
10^ atoms. The six images show the \mp) — component, and increasing magnetic field 
gradient is equivalent to increasing synthetic gauge potential. Reprinted with permission 
of I. B. Spielman. 

ence [50] reviews various proposals for such synthetic gauge fields. Recently, 
Bloch's group in Munich has studied a two-dimensional ^^Rb BEC in a square 
lattice, producing strong effective magnetic fields with staggered (alternating 
striped) fiux of order 1/2 fiux quantum per plaquette [57|. This achievement 
is a significant step toward realizing highly correlated states that are analogs 



of fractional quantum-Hall states [51 



The situation becomes much more interesting if two or more degenerate 
states contribute simultaneously to the Berry's phase [50|, Sec. III]. Instead of 
the previous example of two well-separated states \g) and |e), consider a set of 
q degenerate ground states \Xj)- If the system starts in a linear combination 
of these states, then the adiabatic motion of the particle will remain within 
this degenerate set, leading to the following matrix generalization of the 
synthetic gauge field in Eq. (1^ : 



Ajk = ih{xj{r)\Vxk{r)), 



(68) 



where the effective vector potential is now a. q x q matrix in the internal 
space, as well as a vector in three-dimensional coordinate space. In general. 
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these matrices do not commute, which leads to non-Abehan gauge fields 58 
that are qualitatively different from the Abelian case with q = 1. 

One unusual aspect of these non- Abelian matrix gauge fields appears in 
the effective magnetic field B. By definition, the velocity is the combination 
V = {p — A)/M. Its Heisenberg equation of motion ihi) = [f, |Mv^] leads 
to the expected form Mv = ^{v x B — B x v), where the symmetrized 
structure arises because v and B need not commute. Here, the effective 
magnetic field B is given by the generalized curvature 

1 i 

Bi = -eijkFjk, where Fjk = djAk - dkAj - -[Aj, Ak]. (69) 

Even if these gauge fields are spatial constants, the effective magnetic field 
does not, in general, vanish because of the term ^eijk[Aj,Ak\ = {A x A)i. 
This situation is dramatically different from that of the more familiar Abelian 
gauge field, where the magnetic field vanishes for a spatially constant A. 



Reference [59| analyzes a particular generalization of the experiments by 
the Spielman group involving two degenerate ground states obtained by bal- 
ancing the linear and quadratic Zeeman effect along with the applied field 
gradient. 

Another remarkable feature of these non- Abelian gauge fields is the pos- 
sibility of generating an effective spin-orbit coupling for trapped neutral 



atoms [59|, 160|, |61|. In certain situations, the matrix vector potential A 
can lie in the xy plane and reduce to a linear combination of the 2x2 
Pauli matrices ax and ay. Depending on the vector coefficients, the term 
{p ■ A + A ■ p)/2M in the Hamiltonian can produce a spin-orbit coupling 
involving products of momentum components and Pauli matrices. Simi- 
lar coupling terms p x cr ■ z = Px<^y — Pya^ and Pyay — PxCx are familiar 
in the condensed-matter fields of semiconductor spintronics and topological 
insulators, where they are known as Rashba and Dresselhaus spin-orbit cou- 



pling [62[. In this way, it may be possible to use cold atoms to explore some 
of the properties of these rather complicated semiconductors, although cur- 
rent experiments include only a single Pauli spin matrix and thus provide 
only an Abelian coupling |6l|- There are also proposals for 3x3 spin-1 gen- 
eralizations of the Pauli matrices for spin-i. Reference [50, Sec. III.D and 
Sec. V] summarizes this fascinating field and its connection with the rapidly 
emerging field of spintronics. 
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